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§ ! Abstract 

[ A one-parameter family of time-symmetric initial data for the radial infall of a particle into a 

^2^ ■ Schwarzschild black hole is constructed within the framework of black-hole perturbation theory. 

The parameter measures the amount of gravitational radiation present on the initial spacelike 
surface. These initial data sets are then evolved by integrating the Zerilli-Moncrief wave equation 

m 

K* ' in the presence of the particle. Numerical results for the gravitational waveforms and their power 

^^ ■ spectra are presented; we show that the choice of initial data strongly influences the waveforms, 



o 



both in their shapes and their frequency content. We also calculate the total energy radiated 
by the particle-black-hole system, as a function of the initial separation between the particle and 



O ■ the black hole, and as a function of the choice of initial data. Our results confirm that for large 

Vh . initial separations, a conformally-flat initial three-geometry minimizes the initial gravitational- wave 

content, so that the total energy radiated is also minimized. For small initial separations, however, 

/\^ • we show that the conformally-flat solution no longer minimizes the energy radiated. 

PACS numbers: 04.70.Bw, 04.30.Db, 04.30.Nk 



I. INTRODUCTION 

Due to the advent of laser-interferometric gravitational-wave detectors, there is a strong 
possibihty that gravitational waves from black-hole collisions will be detected in the near 
future. Black-hole mergers are good candidates for detection because of the strength of the 
gravitational waves they emit, and of the relatively high expected event rate[l[]. Detection 
of these collisions and successful extraction of the black-hole parameters require a detailed 
theoretical understanding of the gravitational waves emitted during the collision. The lack 
of solutions to the dynamical two-body problem in general relativity renewed interest in 
simulating black-hole collisions by numerically integrating the full non-linear Einstein field 
equations[|, |^, §, ^ ||, ^ ||, H, |^. The challenge is enormous. For a recent review of the 
field of numerical relativity, see ||TT |. 



From a numerical point of view, long-term simulations of black-hole collisions are limited 
by the available memory^ and instabilities associated with the numerical implementation 
of the full non-linear equations^, H, |l^. Because of these difficulties, numerical relativity 
is not yet at the stage where it can simulate black- hole collisions for a very long timeP|. 
Most simulations are started at a late stage of the collision, when the two black holes are 
separated by a distance of just a few black-hole masses. These small initial separations 
imply that non-trivial initial values must be provided for the gravitational field in order to 
start the numerical evolution. 

The initial- value problem consists of finding an initial three- metric 7jj, and an initial ex- 
trinsic curvature Kij, that encode all of the physical information about the system. The solu- 
tion should contain information about the state of motion of the black holes, and information 
about the gravitational-wave content of the initial three-surface. The initial gravitational- 
wave content is important because it contains information about the motion of the black 
holes prior to the beginning of the numerical evolution. (In the past of the initial hyper- 
surface, the two black holes emitted gravitational waves that must be accounted for in the 



initial data.) As we will see in Sec. p| , the initial gravitational- wave content can play an 
important role in the modeling of the gravitational waveforms. 

For the case of colliding black holes, various methods have been developed to find initial 
data that satisfy the Hamiltonian and momentum constraints of general relativity, e.g., the 
apparent-horizon method[0, conformal-imaging method IT^, and puncture method |[T5[ . 



The starting point in these methods is a maximally embedded initial spacelike slice {K^^ = 
0). On this slice, the geometry is described by a conformally-flat initial metric, jij = ^I/^7ij, 
and a rescaled traceless extrinsic curvature Kij = ^~^Kij^ where 7jj is the metric of three- 
dimensional fiat space, \E' the conformal factor, and Kij the traceless conformal extrinsic 
curvature. With these choices, the momentum constraints in vacuum can be recast in terms 
of quantities defined in the fiat space: V^Kij = 0, where V* is the fiat-space covariant 
derivative. Bowen and Yorlc obtained solutions to this equation |T^ in terms of a vector field 



V^ that contains information about the spins and linear momenta of the holes on the initial 



surface, i.e. K"-^ = 'V'V^ + VW^ - 2/37*-'V V^; explicit forms for V"- can be found in p7 
The three techniques mentioned above rely on the Bowen- Yorlc solutions to construct K^j, 
but because of the different choices of topology for the initial hypersurface, the Kij^s they 
obtain are different. But since the differences in initial topology are mostly hidden behind 
the event horizons of the blaclc holes, they cannot affect significantly the resulting dynamics. 
To reflect this, we will generically refer to a Kij based on the Bowen- York: solutions as a 
"longitudinal extrinsic curvature" . 

Thus, solutions to the initial-value problem can be generated, but the real problem is 
to construct physically suitable solutions. This is a difficult problem because of the non- 
linearities inherent to the theory. It is very hard in practice to identify which components 
of jij and Kij are to be constrained, and which are to be associated with dynamical and 
gauge degrees of freedom. It is now widely accepted that choosing a conformally-fiat ^ij and 
a longitudinal Kij are unlikely to yield physically relevant initial data sets. 

For example, the assumption that the initial metric of a black-hole binary is conformally 
fiat is likely to be inappropriate, because the metric is known not to be conformally fiat at 



the second post-Newtonian order |jT^. As well, any physical black hole is likely to be rotating 



and, in this case, the assumption of conformal flatness yields poor initial data. The reason 



for this is well known: Garat and Price [|lg] have shown that there is no spatial slicing of 
the Kerr solution that is both axisymmetric and conformally fiat, and reduces smoothly to 
the Schwarzschild solution in the no-rotation limit. The techniques mentioned above can 
nevertheless provide solutions to the initial-value problem that represent a rotating black 
hole, but the solutions cannot correspond to a stationary Kerr black hole; there must be 
some gravitational radiation on the initial slice. 

Even in the case of a perturbed Schwarzschild black hole, where conformally-fiat slices 



can be constructed, the solutions obtained with a longitudinal extrinsic curvature are not 



necessarily adequate. Lousto and Price [|T^] have shown that for the head-on collision of 
two nonrotating black holes in which one of the holes is much less massive than the other, 
conformally-flat 7jj and longitudinal Kij data (CFL data) do not reproduce the numerical 
results|T^. In their analysis, they imposed the CFL data at a time to and evolved it forward 
in time. They then looked at the conditions at a later time ti, and found that the extrinsic 
curvature extracted from the numerical data matched poorly with the extrinsic curvature 
obtained at t = ti from CFL data. Instead, they found that the extrinsic curvature was 
better represented by postulating a convective time derivative [pO|] , which means essentially 
that the time derivative of the metric is proportional to the four-velocity of the small black 
hole. 

The various problematic issues associated with the initial-value problem can be better 
understood by having recourse to approximate methods [pl|, p2| . For large initial separations 
and slow-motion processes, post-Newtonian (FN) theory is useful, especially when the FN 
metric of two point masses is matched, in a buffer region, with the metric of two distorted 



Schwarzschild black holes 1 23]. For small mass ratios, black-hole perturbation theory can be 
very useful p^, [25| , |26| , |27[] , since in this case there is no restriction on the velocity of the 
small body. 

This paper is the first of a series in which, following Lousto and Price [P^ pU| , pB[ , the 
initial-value problem will be systematically explored in black-hole perturbation theory. Our 
goal is to explore the space of solutions to better identify physically realistic situations 
(within the domain of validity). This, we hope, will help us understand the essential features 
that should be possessed by a solution to the initial-value problem. In this paper we look at 
the radial infall of a particle starting from rest at a radius Tq outside a Schwarzschild black 
hole, and we construct a one-parameter family of solutions to the initial-value problem. 
The parameter, a, measures the amount of radiation present on the initial hypersurface. 
For a = 1 the initial metric is conformally flat, a choice of initial data that was previously 
considered by Lousto and Price p8[. Starting from our family of initial data, we calculate the 



time evolution of the gravitational perturbations produced by the infalling particle. We then 
calculate the gravitational waveforms, their power spectra, and the total energy radiated. 
In this way, we are able to explore the limitations of the assumption of conformal flatness, 
and determine when the contamination of the gravitational waveforms, due to unphysical 



radiation on the initial three-surface, becomes important. 

The paper is organized as follows. In Sec. |I| we give a short overview of black-hole pertur- 
bation theory, appropriately restricted to the specific case of radial infall into a Schwarzschild 
black hole. The evolution of the metric perturbations, governed by an inhomogeneous wave 
equation, is described in Sec. p.1 A|. Our one-parameter family of solutions to the initial- value 



problem is constructed in Sec. [lllj| . Section |in| is divided into four subsections. In Sec. [ill A 



we explain the numerical method used to integrate the wave equation in the presence of the 



particle. In Sec. |IIIB| we present the gravitational waveforms extracted from the simula- 

we 



tions; we discuss the influence of the choice of initial data on the waveforms. In Sec. [IIIC 



calculate the power spectrum of the waveforms for selected values of a and initial separation 
Tq. In Sec. [1IID| we compute the total energy radiated in the Z = 2, 3, and 4 multipole 



moments of the gravitational radiation field. In Sec. |^ we summarize our findings and offer 
some conclusions. Finally, the Appendix contains a discussion of the convergence properties 
of our numerical code. 

Throughout the paper we use MTW conventions ||29|| for the metric signature and cur- 
vature tensors. Greek indices run from to 3, and Latin indices run over the t and r 
Schwarzschild coordinates. Geometrized units are used: G = c= 1. 

II. PERTURBATION THEORY 



A. Zerilli-Moncrief Wave Equation 

We consider the radial infall of a particle of mass /i, starting from rest at a radius Tq, into 
a Schwarzschild black hole. The particle creates a small perturbation in the gravitational 
field of the black hole and the total metric is expressed as g^^y = g^^ + /i^j,, where g^^ is the 
Schwarzschild metric and h^^, is the perturbation; this can be decomposed into odd-parity 
(axial) and even-parity (polar) modes, each carrying spherical-harmonic indices / and rrapO 



Because of the azimuthal symmetry of the problem, only the polar perturbations are excited 
by the infalling particle. Our calculations are performed in the Regge- Wheeler gauge [pT| in 
which the perturbation — for a specific /-pole and with m = Q — is written as 



h'^Jx^dx'' = Y'^{e) [(1 - 2M/r)W^{r, t)dr + 2W^{r, t)dtdr 



{21) 



where H^fj^, iJ{°, H!f, and K^*^ are the multipole moments of the metric perturbations. Below 
we will omit the superscripts "/O" where there is no risk of confusion. 

Inserting the perturbed metric into the Einstein field equations yields a set of seven 
coupled equations for the four metric perturbations; three of these come from the Hamil- 
tonian and momentum constraints and the other four are evolution equations. (The three 
remaining field equations describe the axial perturbations of the Schwarzschild black hole.) 
The equations are decoupled when written in terms of the gauge- invariant Zerilli-Moncrief 



function 32 



^('^'^) = xTir^'''^^^x[^'^'^'^^~'^¥^^''''^^ 



(2.2) 



where A = (/ + 2)(/ — 1)/2, / = 1 — 2M/r, and A = A + 3M/r. This normalization for ip differs 
from the normalization adopted in [0, ^, |3^. A relation between various normalizations 



used in the literature can be found in Lousto and Price|2^. Our normalization is chosen to 
agree with theirs. 

The evolution of ip{r,t) is governed by a single inhomogeneous wave equation, 



Qj2 Qj.* 2 



V{r) 



^ir,t)=Sir,t), 



(2.3) 



where r* = r + 2M ln(r/2M — 1) is the tortoise radius 

2/ 



V{r) 



r2A2 



2,, , 9M2/. 2M 



A^(A + 1) + 



(2.4) 



is the Zerilli-Moncrief potential, and 

2 



Ar 



A(A - ly + (4A - 9)Mr + ISM^Jg" + r(A - f)Q 



(2.5) 



is the source term appropriate for radial infall. The Q^^'s (a, & = t,r) are constructed from 
the particle's stress-energy tensor: 



Q^\r, t) = 87T f T''\r, t, 6, y^)Y*{e, (^)dfi. 



(2.6) 



where T"^ = jj, j u°'u^5^{x'^ — Xp(r)) {—g)^ dr, r is proper time along the particle's world 
line Xp^r), u" = [E/f, —{E"^ — fY^'^, 0,0] is the four-velocity of the particle, and E is the 
conserved energy per unit mass. We assume, without loss of generality, that the motion 



proceeds along the negative z direction, so that 6' = on the world line. Evaluating the 
integrals in Eq. ( p.6| ) and substituting the results into Eq. ( |2.5| ), we obtain^ 



S(T,t) 



,,.^^^^l^t!lU(r-rM 



l{l 
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Ar2/ 



5{r-T,m"\ 



where rp{t) describes the radial component of the particle's motion, a " ' " denotes a deriva- 
tive with respect to r, and we have used y'°(0) = J {'2,1 + l)/(47r), which is appropriate for 



motion along the 2;-axis. For the radial infall of a particle starting from rest, Tp(t) is given 
implicitly by 
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AM 



(2.8) 



where Tq = rp(t = 0) is the initial position of the particle, and E = Jl — 2M/ro. In Sec. |T| 
we numerically integrate Eq. ( |2.3| ) with the source term of Eq. ( |2.7| ). The integration of 
Eq. (|2.3|) requires the specification of initial data for the Zerilli-Moncrief function: both 



■ip(r,t = 0) and ^■ip{r,t = 0) must be known. We shall examine this issue in the next 
subsection. 

The Zerilli-Moncrief function ?/'(r, t) is directly related to the gravitational waves received 
at infinity^, and knowledge of ip{r,t) allows us to calculate the rate at which gravitational 
waves carry energy to infinity. For a specific /—pole, the luminosity (the gravitational- wave 
flux integrated over a two-sphere) is given by 

2 

I , (2.9) 



1 {l + 2)\f^^^ 



±E - 

du ' 647r(/-2)! \du 



where u = t — r* is retarded time. The total energy radiated is obtained by integration: 



1 (/ + 2)! foo ( d ' 
^' = 64^(7^2)! i-oo^H^^ 



(2.10) 



^ The source term given here agrees with the source term given in ||l9|, ^ if we substitute "(1 — 
2M/rp{t))d'{r - rp{t)) - 2M/r'^S{r - rp(i))" for "(1 - 2M/r)S'{r - rp(t))" in our source term. 

^ With polarization axes oriented along the 9 and ip directions, wc have that the "plus" polarization of 
the gravitational field is given by h+ — ^\/{l + 2)(/ + 1)1(1 — l)^p{u) -2Y^^{9, ip)/r, where -2^'" are spin- 
weighted spherical harmonics |34|. The "cross" polarization, in this orientation, is zero. 



B. Initial Data 

The initial-value problem of general relativity is solved once the metric 7jj and the ex- 
trinsic curvature Kij are specified on a three-dimensional spacelike hypersurface. These 
quantities are not freely specifiable, because they must satisfy the Hamiltonian and mo- 
mentum constraints of general relativity. Because the constraints are non-linear, finding 
solutions to the initial-value problem is a non-trivial task, and there is no unique way of 
doing this. Another essential difficulty is to find solutions that represent a physical situation. 
In the case of black-hole collisions, 7jj and Kij must represent two moving black holes with 
some initial amount of gravitational waves, as was explained in the Introduction. 

Two much-studied solutions to the initial-value problem are the Misner solution |]35[ 



(which is generalized by the conformal-imaging method to situations where the initial hy- 
persurface is not a moment of time-symmetry) and the Brill-Lindquist solution [^] (which 
is generalized by the puncture method). They represent the spatial metric of two black 
holes momentarily at rest, and about to undergo a head-on collision. Because the initial 
hypersurface is a moment of time symmetry, K^j = and the initial- value problem reduces 
to finding 7,^ (or the conformal factor). These solutions and their generalizations were used 
as initial data for the study of the head-on collision of two equal-mass black holes in full 
numerical relativity 0, |^, |26|] . 

For the problem studied in this paper, namely the radial infall of a small mass /i into 
a Schwarzschild black hole, we will concentrate on Brill-Lindquist (BL) type initial data. 
(Our initial data sets are not taken to be inversion symmetric with respect to the throat of 
the black hole.) The Brill-Lindquist solution can be expanded in powers of fi/M to yield 
an approximate solution to the initial-value problem appropriate for perturbation theory; in 
this case Hi = and 



v/4vr/(2/ 



H2 = K = 2a- V ; , ' , < , (2.11) 

^(l + M/2fo)(l + M/2f)f^+i' ^ ^ 

where f = r(l + ^7)^/4 is the isotropic radius, / = 1 — 2M/r, and f^ (f>) is the smaller 
(greater) of f and fg. This solution was first examined by Lousto and Pricep8|]. Here we 



wish to construct a more general solution to the initial-value problem in perturbation theory. 

Our initial situation is that of a small mass fi momentarily at rest at a radius Tq outside a 

Schwarzschild black hole; this is a moment of time symmetry, which implies that K^j = 0, so 



that the momentum constraints are automatically satisfied. To complete the solution to the 
initial-value problem, we write the Hamiltonian constraint ^R = Svrp {^R is the Ricci scalar 
on the initial hypersurface, and p is the matter density) in terms of the metric perturbations 
in the Regge- Wheeler gauge (Eq. ( |2.1D ): 

- r^f^K - (3r - 5M)^K + rf-^H2 + [H^ - K) 



dr^ dr dr 



+-^^{H, + K) = r'fQ^' = 8n^^pE5{r - r„). (2.12) 

This is a single equation for two unknowns: the metric functions K{r) and H2{r). One of 
the functions can be specified freely, and we do this by postulating the linear relationship 
H2 = <yK, where a is a constant. This is motivated by the fact that the relation H2 = K 
produces a metric that is conformally fiat: 

= (1 + Kio{r)M0)) (1 + — j (df2 + fMQ') . (2.13) 

Making the choice H2 = K turns Eq. ( p.l2| ) into a hypergeometric equation for K{r), and 
as we shall show below, its solution is given by Eq. (|2.11 ) above. The more general relation 



H2 = aK also turns Eq. ( p.l2|) into a hypergeometric equation, for which we will be able 
to find solutions that generalize Eq. ( p.llj ). Our solutions are parameterized by a, and the 
procedure outlined here produces a one-parameter family of initial-data sets. Setting a = 1 



reproduces the BL solution of Eq. (|2.11 



Postulating the relation H2 = aK at the initial moment t = 0, Eq. ( ^.12] ) becomes 

z = r/2M) 

d' 



'^'-^h.^- 



;K + (5/2 - « - (3 - a)z) ^K - ({1 -a)- ^^^^(1 + a)) K = 0,(2.14) 

on either side of the particle's initial position (r 7^ Vo). This is the hypergeometric equation. 
For r < To we must choose a solution that is regular at r = 2M, but that is allowed to 
diverge at r = 00. We denote this solution by K^{r). For r > Tq we must choose a solution 
that is regular at r = 00, but that may diverge at r = 2M. This solution is denoted by 
Ky{r). In terms of these, the solution for K{r) takes the form^ 

K{r, 0) = C(r„)ir<(r<)ir>(r>), (2.15) 



■^ In principle, we are free to add any multiple of Ky (r) to our particular solution. This represents a different 
choice of initial three-geometry ]19| , and we do not consider this possibility here. 



where r< (r>) is the smaller (greater) of r and Vo, 

-b 



K<{r) 
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2MJ i^(M-c+l;l-a + 6;^) 



2Af 
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2M 
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(2.16) 



F(, ; ; ) is the hypergeometric function, 



a 1 

a = 1 h - 

2 2 

a 1 



2 2 



0^ + 2(1 + «)/(/ + !) 



a2 + 2(l + a)/(/ + l) 



1/2 



1/2 



c = a, 

2 



(2.17) 



and C{ro) is a constant ensuring that ^K has the correct discontinuity at r = Tq. Integrating 



across the (5-function in Eq. ( p.l2| ), using H2 = aK and Eq. ( p.l5|) for K^ C{ro) is easily 
found to be 



r( \ 8 ^ 2/ + 1 1 r 
C[ro) = 8tc-~ 



E\ 47r rl 



K^{r,)K^{r,)~K^{r,)Ky{r,) 



47n 



121 + 1 fx (To \^'''T{ro, 



4:71 M \2M 



E 



(2.18) 



where 



Tir, 



aF{a + l,a-c+l;l-6 + a; 2M/ro)F{b, b-c+l;l-a + b; 2M/ro) 



bF{a, a-c+l;l-b + a; 2M/ro)F{b + 1, 6 - c + 1; 1 - a + 6; 2M/ro) 



(2.19) 



The second equality in Eq. (|2.18| ) follows from Eq. (15.2.2) of |^ together with Eq. (|2.16|) 
to evaluate K_^{ro) and K-^ijo). 

The initial value of the Zerilli-Moncrief function can now be obtained by inserting K{r) 
and H2 = aK(r) in Eq. (|2.2|). This, together with ^'?/'(r, 0) = 0, completely specifies our 
family of initial-data sets for the radial infall of a particle starting from rest. In Sec. |T| we 
numerically integrate the Zerilli-Moncrief equation for this family of initial data. 

The conformally-flat solution can easily be recovered as a special case of the initial-data 
sets presented above. For a conformally-flat initial three-geometry (a = 1), the parameters 
appearing in the hypergeometric functions (Eq. ( ^.17| )) take the values a = I + 1, b = —I, 



10 



and c = |, and Eq. ( |2.16|) becomes 



K. 



:W = (^)'FH,-/-l/2;-2/;2M/r) 

/ r \ I / I — \ 21+1 



K>{r) 



2Mj l + M/2f 



(2.20) 



2M 



F(/ + l,/ + l/2;2(/ + l);2M/r) 



-,21+1 I "^ 



.2M 
/2My+^ 1 



C+l) / ^\ -(2/+1) 



M/2f^ 



(2.21) 



where we have used Eq. (15.1.13) of ||3^. With these results, it is a trivial matter to show 
that 

2M fl 



K4r^)Ky{ry) 



(l + M/2fo)(l + M/2f)f 



i+i- 



(2.22) 



With the conformally-flat values of a, b, and c, r(ro), given by Eq. ( p.l9| ), becomes 



rfr. 



(/ + 1)F (/ + 2, / + 1/2; 2(/ + 1); 2M/ro) F (-/, -/ - 1/2; -2/; 2M/r,) 



n -1 



+/F (-/ + 1, -/ - 1/2; -21; 2M/ro) F{l + l,l + 1/2; 2{l + 1); 2M/ro) 



E 



21 + 1 

E 
2/ + 1' 



F(-Z, -Z - 1/2; -2/; 2M/ro)F{l + 1, / + 1/2; 2(/ + 1); 2M/ro) 



(2.23) 



where the second equality follows from the second identity of Eq. (15.1.13) of |^ applied 
to F{1 + 2, ; ; ) and F{1 + 1, ; ; ), and the third equality from the first identity of Eq. (15.1.13) 
of 1^ applied to both hypergeometric functions. Substituting this results into Eq. (|2.18|) 
yields 



C{ro) 



Att ji 



V 2/ + 1 M 
Finally, inserting Eq. ( 2.22 ) and Eq. ( 2.24 ) into Eq. ( ^.15| ) gives 



(2.24) 



K{r) = 2^1- 



477/(2/ +1) 



r^ 



l + M/2fo)(l + M/2f)f> 
11 



«+i' 



(2.25) 
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FIG. 1: Initial values of ip{r,t) for / = 2; the particle is initially located at ro/2M = 40 (left) and 
ro/2M = 2 (right). In both cases V'('^;0) is peaked and discontinuous at the particle's location, 
r = Vq. Decreasing (increasing) the value of a increases (decreases) the amplitude and width of the 
initial pulse. These properties can be associated with the amount of gravitational radiation present 
on the initial hypersurface. For fixed q, the amplitude of the peak decreases with decreasing Tq. 



which is indeed the conformally-flat solution of Eq. ( |2.11| ) . 

In order for the perturbed metric to be real (as opposed to complex), the solutions to 
Eq. ( |2.14| ) must be real functions. Thus, a must be such that the parameters a and 6, as 
listed in Eqs. ( |2.171 ), are real. It is easy to verify that the quantity appearing under the 
square root will be positive if a > a+, or a < a_, where 



«± = -/(/ + 1) ± J {I + 2)(/ + !)/(/ - 1). 



(2.26) 



It can be verified that a+ varies between —1.10102 (when / = 2) and —1 (for / -^ oo). On 
the other hand, a_ monotically decreases from —10.8990 (when / = 2) as / increases. Hence, 
the metric functions will be real for all values of / if a is restricted by a > —1. 

Figure |l| displays the Zerilli-Moncrief function ?/'(r, 0) for the cases VojIM = 40 and 
ro/2M = 2. The figure shows that ^(r, 0) is peaked at the particle's location, and its 
amplitude and width change with the value of a; the function is discontinuous at r = 



To because of the discontinuity in the term dK/dr appearing in Eq. (2^). We see that 
the amplitude of the peak increases with decreasing values of a, and this effect is more 
pronounced as Tq decreases. Intuitively, we associate a smaller (larger) amplitude with a 
smaller (larger) amount of gravitational radiation initially present in the spacetime. Thus, 



12 



to initial data with a small (large) value of a we associate a large (small) amount of initial 
radiation. 

III. NUMERICAL RESULTS 
A. Numerical Method 

We now describe the numerical method used to integrate Eq. (|2.3|) , starting with the 
initial data constructed in the preceding section. We recall that the source term is given 
by Eq. ( |2.7| ); special care must be given to the fact that S{r,t) is singular at r = rp{t), the 
trajectory of the infalling particle. A numerical method was presented by Lousto and Price 
||T9| and we employ it here with a slight modification, which we discuss below. 

The numerical domain is a staggered grid in r* and t, with a stepsize A = Ar*/4M = 
At/2M, over the region of spacetime bounded by the spacelike hypersurface t = (the 
moment of time symmetry), a null hypersurface approximating the event horizon [located 
at M = T - r*{T), where T is given by Eq. (^ with rp{T)/2M = 1.0001] ^, and a null 
hypersurface approximating future null- infinity [located at v/2M = (t + r*)/2M ^ 1500]. 
The Zerilli-Moncrief function is extracted from the numerical data on this last hypersurface, 
i.e. ip{r,t) is evaluated at v/2M ^ 1500 and expressed as a function oi u = t ~ r*. The 
results presented in Sec. [1IIB| were obtained by setting A = 0.01 in the evolution scheme 



presented below. 

The finite-difference method must take into account the source term, which is non zero 
(and singular) on the world-line of the infalling particle. The numerical grid is divided 
into cells of area 4A^. (Instead of using the physical area 2A^, we found it convenient, 
following Lousto and Price||T^, to use A = J dudv = 4A^.) The cells can be separated into 
two groups. The first group corresponds to cells in which r 7^ rp{t) everywhere; these cells 
are never traversed by the particle. The second group is such that r = rp{t) somewhere 
within the cell; these cells are traversed by the infalling particle. The two types of cells are 
displayed in Fig. ^; in the diagram, the particle enters the cell on the right of r^, and leaves 



^ This choice is made to ensure that the particle's contribution to the radiation is (almost) zero at the end 
of the numerical integration. This means that the value of u associated with the event horizon in our 
numerical grid changes every time we change Tq. 
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(r*,t^+2A) 



(r*,t^+2A) 




(Rj, Tj) 



(r* + A, t + A) 
b D 



(r* + A, t + A) 
b b 




(Rb- V 



FIG. 2: The left cell is never crossed by the particle's trajectory. In such a cell the evolution of 
■0(r, t) is unaffected by the source term in Eq. ( |2.3| ), and Eq. (|3.lD is used to evolve ■0(i,r) forward 
in time. The right cell is traversed by the particle's world line. The areas shown in the cell are 
used in the numerical algorithm that incorporates the effect of the source term in Eq. (p.3|). The 
wave function is evolved according to Eq. (^). The labels (i?^,Tfc) and {Rl,Tt) refer to the points 
{r*,t) at which the particle enters and leaves the cell, respectively. 



on the left. 

The evolution of ip{r* ^t) across cells belonging to the first group is not affected by the 
source term. For these cells, we use a standard scheme, accurate to (9(A^), for the homoge- 
neous wave equation in the presence of a potential: 

A^ 



V'(r,*,tfe + 2A) = -V;(r*,tt,)+ ^(r* + A, t^ + A) + ^(r* - A, t^ + A) 



K 



(3.1) 



where (r^, t^) designates the bottom corner of the cell (see Fig. ^, and VJ, = V{rl). 

The evolution of ip{r,t) across cells belonging to the second group is affected by the 
singular source term. To obtain the evolution of if) across these cells, we closely follow 
Lousto and Price|]19[, carefully keeping terms up to order (9(A^). Integrating Eq. (|2.3|) over 
a grid cell, term by term, we get: 



dA -— r + 



Qf2 Qj.*2 



-^(r* - A, tfc + A) - ij{rl + A,h + A)] 



(3.2) 



dAVir*)i; = Virl) {AMr; - ^,h + A) + A2ijirlh) 

+ Asij{rl h + 2A) + A^^{r; + A,h + A)) + 0{A'), (3.3) 
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where dA = dudv, and the Aj's are areas dividing the cell into four non-equal parts (as 
shown in Fig. |^) . The integration of the source term over the cell eliminates the (5-function in 
Eq. ( p.7| ) and the term involving the derivative of the ^-function is evaluated by integrating 
by parts. The result is 

'6M 



iAs = -.r-itm. 



±K 



n 

i fin) 
[Hn] 



rp{t) 



'l-E') + \{\ + l) 
fiTt) 



3M2 ^ M 

+ 4A- 



rAty 



rp{t) 



[lTr;(T,ri + ^[l±r;(T,ri 



(3.4) 



where k is defined by 1{1+1)Ek = 167r^^(2/ + l)/(47r), r; = -^JE^ - f/E, f{t) = f{rp{t)), 
A(t) = A(rp(t)), Th is the time at which the particle enters the cell, and Tf the time at which 
it leaves the cell. In the previous expression, the upper (lower) sign for the first boundary 
term (a function of Tb) applies when the particle enters the cell on the right (left) of r^. 
Similarly, the upper (lower) sign for the second boundary term (a function of Tt) applies 
when the particle leaves the cell on the right (left) of r^. 

Substituting the previous results in Eq. (|2.3|) and solving for ip^r^, U + 2A), we obtain 



^(r*,4 + 2A) 



1 + ^(^2-^3) 



+V'(r* + A,tfc + A) 



1-^(^1 + ^3) 
Vb 

1-^(^4 + ^3) 



+V'(r*-A,4 + A) 

A-^ j j dAS{r,t) 



-\- 



4 



(3.5) 



Notice that we have kept the factor of (1 — ^^3) in front of the source term. This factor 
does not appear in the scheme devised by Lousto and Price 0], but is required to obtain 
quadratic convergence: the boundary terms appearing in Eq. (|3.4| ) are 0{\) in A and ^A^ 
represents a correction of order (9(A^) to their contribution, which we cannot neglect if we 
want quadratic convergence. 

The previous equations cannot be used to evolve if) from the initial time t = to the 
next time t = A. To get '?/'(r^,A) we need ip{r^ — A,0), '?/'(r^ + A,0), and ■ip{rl,—A). 
By construction, tp{rl, —A) is not known, but we can use the simple relation ^(r^, —A) = 
ip{rl, A) + 0{A^) which holds by virtue of the fact that the data at t = is time-symmetric: 
dip/dt = 0. This relation is valid so long as r*{t) 7^ r^ everywhere in the cell. We chose 
rl — A < r*(0) < rl which, for radial infall, guarantees that the relation holds. 
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The gravitational waveform produced by the infalhng particle is given directly by '?/'(«), 
which is extracted from the numerical data in the way described previously. Other quantities 
of interest are the energy radiated and its power spectrum. We calculate these by first taking 
a fast Fourier transform of the wave function ^(n). This is given by 

^{uj) = / dcue^'^^VH, (3.6) 



where u is the frequency. From Eq. (|2.9| ) and Parceval's theorem, we find that the power 
spectrum is 

The total energy radiated in each multipole moment is calculated by performing a Romberg 
integration over all frequencies. 

B. Waveforms 

The radiation arriving at the null boundary of our domain of integration contains ra- 
diation from the initial data, radiation emitted by the particle as it falls toward the black 
hole, and radiation corresponding to the black-hole's response to the perturbation. Because 
the particle is at rest at t = 0, it initially produces very little radiation. The early part 
of the waveform is therefore dominated by radiation contained in the initial data. As time 
proceeds, the radiation produced by the particle becomes noticeable, and the dynamics asso- 
ciated with the perturbation of the event horizon starts to play a role. The radiative process 
can therefore be separated into three stages. The first stage is associated with the initial 
data, the second with the infalling particle, and the third with the event-horizon dynamics. 
We now describe these stages in detail. 

The initial data is time-symmetric, and it therefore consists of two radiation pulses: one 
pulse is outgoing, and the other is incoming. The outgoing pulse proceeds to infinity, with 
some backscattering [which is small, unless the pulse originates in the strong-field region of 
the spacetime (r* < 0)]. On the other hand, the ingoing pulse moves toward the black hole 
and is backscattered by the potential barrier at r* ^ 0. The reflected pulse then proceeds 
to infinity, where it arrives delayed with respect to the original outgoing pulse. This is 
the first stage of the radiative process, and it is directly associated with the initial data. 
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FIG. 3: The / = 2, 3, and 4 modes of the Zerilh-Moncrief function for infall from ro/2M = 40, for 
a = —0.3, 0, 1, 2, and 5. For each multipole moment, the early-time behaviour is dominated by the 
initial data contribution {u/2M < 350), while the late-time behaviour is dominated by radiation 
emitted by the particle in the strong-field region (350 < u/2M < 400) and by quasi-normal ringing 
of the black hole {u/2M > 400). The early-time portion of the waveforms (first stage) depends 
strongly on the parameter a, which labels the choice of initial data. In contrast, the late-time 
portion (second and third stages) are completely insensitive to the choice of initial data. 

The second stage is associated with the motion of the particle. As the particle accelerates 
toward the black hole, it produces radiation which propagates to infinity. This happens 
either by direct propagation or by backscattering from the potential barrier. The third and 
final stage of the radiative process is the response of the black hole to the perturbation. As 
the particle approaches the black hole, it tidally deforms the event horizon, which becomes 
dynamical. The radiation produced in this process interacts strongly with the potential 
barrier outside the black hole, and the result, after transmission to infinity, is a pattern of 
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damped oscillations. This response of the black hole to the perturbation created by the 



particle is known as quasi- normal ringing 1 38 



As long as Vo ^ 2M, the different stages of the radiative process can be clearly distin- 
guished (initial-data pulses, acceleration radiation, and quasi-normal ringing). But when Tq 
is comparable to 2M, the three epochs become confused, and this gives rise to interfering 
waveforms. For this reason, varying the initial separation between the particle and the black 
hole can have an important effect on the gravitational waveforms. 

Let us describe more fully the first stage of the radiative process. For infall from a large 
distance {vo ^ 2M), the outgoing pulse travels directly to future null-infinity with very 
little backscattering, since the pulse originates in a region where the potential V{r) is weak. 
This is shown in the waveform as a single pulse of radiation at early times {u ~ —r*). The 
ingoing pulse, on the other hand, proceeds toward the black hole, where it is almost entirely 
backscattered by the potential barrier. The reflected pulse then proceeds to infinity and 
this gives rise to a second pulse of radiation at m ~ r*. Varying a changes the amount of 
gravitational radiation initially present at t = 0, and the amplitude of the two pulses depends 
on a. This can be seen in Fig. ^, where the / = 2, 3, and 4 modes of the Zerilli-Moncrief 
function are displayed for infall from ro/2M = 40, for several values of a. For / = 2, the two 
pulses have a minimum amplitude when a = 1, and the amplitudes vary smoothly with a. 

At later times, once the pulses have made their way to infinity, the radiative process 
becomes dominated by the particle's contribution; this is the second stage. At this time, 
the particle has entered the strong-field region of the spacetime, the acceleration is large, 
and it radiates strongly. This stage lasts for a short time, because the particle quickly falls 
into the black hole. The burst of radiation from the particle is then quickly replaced by 
the third stage, quasi-normal ringing. When ro ^ 2M , the acceleration radiation and the 
quasi-normal ringing stages are insensitive to the choice of initial data. This is illustrated 
in Fig. ^. 

The situation changes when Vo is chosen within the strong-field region of the spacetime. 
The outgoing pulse from the initial data still proceeds directly to infinity with little backscat- 
tering, creating a pulse at u ^ — r* (see Fig. ^. However, the ingoing pulse is no longer 
entirely backscattered by the potential barrier; part of the pulse is now transmitted to the 
black hole. The backscattered portion of the ingoing pulse proceeds to infinity where it gen- 
erates a second pulse at u ~ r*, while the transmitted pulse reaches the event horizon. As a 
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FIG. 4: The / = 2, 3, and 4 modes of the Zerilli-Moncrief function for infall from ro/2M = 5, for 
a = 0, 1, and 2. The waveform changes smoothly as a is varied away from 1 (conformal flatness), 
but the change is more dramatic when a is decreased. The three stages of the radiative process now 
overlap and become confused. However, we still witness an early-time sensitivity, and a late-time 
insensitivity, to the choice of initial data. 

result, the event horizon becomes distorted, and starts radiating into quasi- normal modes; 
the second pulse is therefore immediately followed by an epoch of quasi-normal ringing ex- 
cited by the transmitted pulse. The amplitude of the quasi-normal ringing depends on the 
strength of the excitation, and is therefore highly sensitive to the choice of initial data. This 
early epoch of quasi-normal ringing cannot be seen for infall from large Tq, because in that 
case the initial data contains mostly low-frequency gravitational radiation that is almost 
totally reflected by the Zerilli-Moncrief potential ||39||. 

For intermediate values of Tq (4 ^ ro/2M < 10), the particle starts in the strong-field 
region of the spacetime, and its large acceleration causes it to radiate strongly almost im- 
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FIG. 5: The / = 2 and 3 modes of the Zerilli-Moncrief function for a particle faUing from ro/2M = 
2.1, for a = 0, 1, and 2. At intermediate times, the waveform is a superposition of the backscattered 
ingoing pulse from the initial data, and acceleration radiation from the particle. The choice of initial 
data has an influence at late times; it affects the amplitude and the phase of the quasi- normal 
ringing. 

mediately. The second stage, which for infall from large Tq was dominated by acceleration 
radiation, is now the sum of acceleration radiation and quasi-normal ringing excited by the 
initial ingoing pulse, as was discussed above. This can be seen especially clearly for the / = 3 
and 4 modes of the Zerilli-Moncrief function displayed in Fig. ^, for infall from ro/2M = 5. 
The superposition can be seen in the interval 5 < u/2M < 20; it is small for a = 1, but large 
for a 7^ 1. The third stage of the radiative process is the quasi- normal ringing of the black 
hole, excited by the particle as it reaches the event horizon. This phase of quasi-normal 
ringing is to be distinguished from the earlier phase associated with the transmitted pulse; 
this new phase is insensitive to the choice of initial data. 

As To is decreased further (1.3 ^ ro/2M ^ 4), the waveforms become increasingly con- 
fused. In this range of initial separations, the particle radiates strongly immediately, and 
the reflected pulse of ingoing radiation does not reach infinity before the radiation from the 
particle becomes significant. The interference between these two contributions to the wave- 
form is fairly small, because the particle radiates for a short time before passing through 
the event horizon. Figure ^ displays the / = 2 and 3 modes of the Zerilli-Moncrief function 
for infall from ro/2M = 2.1 for a = 0, 1, and 2. The interference between the backscattered 
ingoing pulse and the radiation emitted by the particle is apparent, especially for / = 2 and 
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FIG. 6: The / = 2 and 3 modes of the Zerilli-Moncrief function for a particle faUing from ro/2M = 
1.1, for a = 0, 1, and 2. For this very small initial separation, the particle is quickly absorbed by 
the black hole and cannot radiate much. The initial value of ip{r,t) determines how strongly the 
event horizon is distorted, and it therefore has a major impact on the quasi-normal ringing phase 
of the radiation. 

a = 0. A new phenomenon is observed for these values of Tq: the stage of pure quasi-normal 
ringing is now affected by the choice of initial data. The event horizon is distorted both by 
the transmitted ingoing pulse and the particle, and these factors act at roughly the same 
time. Different choices of initial data will therefore affect differently the amplitude and the 
phase of the quasi-normal ringing. This is displayed in Fig. |^. 

These effects disappear when To is moved past the potential barrier (1 < ro/2M ^ 1.3). In 
such cases, most of the outgoing pulse in the initial data is reflected by the potential barrier, 
and does not register at infinity. Instead, the pulses in the initial data and the radiation 
from the particle work together to distort the event horizon. For these very small Tq, the 
quasi-normal ringing is the only feature that remains in the waveforms. The amplitude and 
the phase of the quasi-normal ringing are very sensitive to the choice of initial data, because 
the information about the initial distortion of the event horizon is entirely encoded in the 
initial data (see Fig. ^ for the case ro/2M = 1.1). The contribution from the particle is 
minimal because it is almost immediately absorbed by the black hole. 
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FIG. 7: Power spectra for the / = 2, 3, and 4 modes of the Zerilh-Moncrief function for ro/2M = 40, 
and a = 0, 1, and 2. Varying a changes the shape of the spectra at fow frequencies, but the effect 
disappears at high frequencies. The low-frequency oscillations are more important for a = 0, 
compared with a = 1 and 2, and a = 2 produces more oscillations than a = 1. The low- frequency 
part of the spectrum is due mostly to the low-frequency gravitational waves contained in the initial 
data. 

C. Power Spectra 

The power spectra, as calculated from Eq. ( p.7|) with a fast Fourier transform algorithm, 
tell a similar story. For a given choice of /, a, and Tq, the total power spectrum is the sum 
of initial-data, particle, and quasi-normal ringing contributions, but it includes interference 
between these contributions. 

The interference is important when the pulses associated with the initial data, the ac- 
celeration radiation, and the quasi-normal ringing contain overlapping frequencies. Typical 
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frequencies are lUq ~ r„ ^ for the initial-data radiation, tUp ~ r^ ^ for the acceleration radi- 
ation, and ujQN ~ M~^ for the quasi-normal ringing. As was mentioned in Sec. [Ill B| , the 



particle emits mostly in the strong- field region of the spacetime {ujp ~ M~^) and conse- 
quently, the interference between acceleration radiation and quasi-normal ringing is always 
present, while interference with the initial-data contribution is important only when Tq is 
not much larger than 2M. In other words, interference effects involving the initial-data 
pulses are important when the particle is released in the strong-field region. This picture is 
somewhat simplistic, but it serves as a useful guide to determine when interference effects 
become important. Our numerical results are consistent with this picture. 

For large initial separations (ro/2M ^ 1), we expect the energy spectrum to be the direct 
sum of the powers in the initial-data and particle contributions, and in the quasi-normal 
ringing, without much interference. This is confirmed in Fig. |^, where the spectra for the 
/ = 2, 3, and 4 modes of the Zerilli-Moncrief function for ToJIM = 40, and a = 0, 1, 2, 
are presented. At low frequencies, the initial data manifests itself as a strong pulse which 
dominates the spectrum. At higher frequencies, acceleration radiation and quasi-normal 
ringing dominate the spectrum; this part of the spectrum is easy to recognize, as it does not 
change when a is varied. As we move away from conformal flatness (a = 1), the influence 
of the initial data spreads into higher frequencies. This can be seen as oscillations in the 
spectrum (cf. the cases a = and 2 in Fig. |^). These oscillations are also present for a = 1, 
but increasing or decreasing the value of a increases their amplitude and the extent by which 
they spill into higher frequencies. In general, a choice of initial data with a < 1, instead 
of a > 1, produces larger oscillations, and the effect extends to higher frequencies. (The 
oscillations are the smallest when a = 1.) 

We have seen in the previous section that for intermediate values of Tq (4 ^ ro/2M < 10), 
the choice of initial data affects the shape of the waveforms up to times where acceleration 
radiation starts to dominate the radiative process. For these initial separations, interference 
effects are important. This is confirmed by our numerical simulations of infall from ro/2M = 
5. The spectra, for this value of Tq, are displayed in Fig. || for / = 2, 3, and 4, and a = 0, 
1, and 2. For 1 = 2, initial data with a > 1 tend to amplify the features apparent for 
a = 1 (see the case a = 2 in Fig. ||), while decreasing the value of a changes the location of 
the maxima; local minima appear in the spectrum where maxima were seen for a = 1 (see 
a = in Fig. ^. For / = 3 and 4 and a > 1 (a < 1), we observe a similar amplification 
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FIG. 8: Power spectra for the / = 2, 3, and 4 modes of the Zerilh-Moncrief function for infah from 
ro/2M = 5, and a = 0, 1, and 2. Interference between the initial-data pulses, the acceleration 
radiation, and the quasi-normal ringing plays a crucial role in determining the shape of the spectra. 
For a = and / = 2, the spectrum has two maxima at 2Mlo = 0.4065 and 2Muj = 0.6455, 
frequencies at which the spectra for a = 1 and 2 have a local minimum. Similarly, the maximum 
at 2Muj = 0.5364 for a = 1 and 2 is replaced by a minimum when a = 0. This is indicative 
of strong interference effects. The labels "xlO" and "xlOO" indicate the amount by which the 
amplitude of these two spectra were multiplied to be presented in the same figure. 

(attenuation) of the features present for a = 1, but the effect is much weaker than for 1 = 2. 
As To is taken closer to the potential barrier (1.3 ^ ro/2M ^ 4), the interference becomes 
less important. In Sec. [111IJ| , we showed that for Tq close to the potential barrier, the particle 
does not radiate strongly before passing through the event horizon, and the interference is 
small because of the small amount of acceleration radiation. In Fig. |^ we present the spectra 
for infall from ro/2M = 2.1, a = 0, 1, and 2, and I = 2, and 3. Although interference effects 
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FIG. 9: Power spectra for the / = 2 and 3 modes of the ZerilU-Moncrief function for infall from 
ro/2M = 2.1, and a = 0, 1, and 2. For / = 2, we also display a = 0.35, to show the smooth transi- 
tion between the single-maximum and two-maxima regimes described in the text. For such a small 
separation, acceleration radiation is small, and the interference effects are not very pronounced. 
They still, however, play an important role in determining the shape of the spectra. The spectrum 
for I = 2 and a = has a shape that indicates strong interference effects between initial-data, 
particle, and quasi-normal ringing contributions: the single maximum at 2Muj = 0.6553, seen for 
a = 1 and 2, is replaced by two maxima at 2Muj = 0.4915 and 2Muj = 0.8697. 

are not as important as for infall from ro/2M = 5, they still play a role in determining the 
shape of the spectrum, as can be seen for the case a = 0. In this case, the spectrum has two 
maxima and a single minimum; this minimum occurs close to the fundamental quasi-normal 
frequency, which suggests a strong destructive interference between initial-data-excited and 
particle-excited quasi-normal ringing. This is different from the cases a = 1 and 2, for 
which the spectra contain a single peak, indicating that most of the energy is radiated into 
quasi-normal modes. 

Interference effects, such as the ones shown in Fig. P, have been observed previously by 



LoustoPDI, but in a different context. Instead of evolving the perturbed Misner solution, 
appropriate for a head-on collision in perturbation theory, Lousto chose to use the full Misner 
solution [^ as initial data for the Zerilli-Moncrief function. He then evolved these initial 
data using Eq. ( p.3|) . In this case, the interference is due to the non-linear nature of the 
initial data he evolved. 

Finally, when Tq is chosen well inside the potential barrier (1 < ro/2M ^ 1.3), the 
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FIG. 10: Power spectra for the / = 2 and 3 modes of the Zerilh-Moncrief function for infall from 
Tq/'IM = 1.1, and a = 0, 1, and 2. For Tq well inside the potential barrier, acceleration radiation 
is negligible since the particle is absorbed by the black hole immediately. The initial data distorts 
the event horizon, which becomes dynamical and starts radiating at its quasi-normal frequencies. 
The amplitude of the quasi-normal ringing is determined by the strength of the tidal distortion 
exerted on the event horizon by the initial data. For the cases displayed, the lowest amplitude is 
obtained for a = 1, while the highest is obtained for a = 0. 

interference effects mentioned previously disappear. For small Tq, the magnitude of the tidal 
distortion applied to the event horizon is affected only by the choice of initial data. Since the 
radiation is then dominated by quasi-normal ringing created by the initial tidal distortion 
of the black hole, the spectra contain a single peak, and its position is independent of a. 
The choice of initial data affects only the amplitude of the quasi-normal ringing, and this 
is reflected in the amplitude of the spectra at the fundamental quasi-normal frequency. For 



infall from ro/2M = 1.1, displayed in Fig. [1^ for / = 2 and 3, and a = 0, 1, and 2, the lowest 
amplitude occurs for a = 1, while the highest occurs for a = 0. 



D. Total Energy Radiated 

In this subsection we calculate the total energy radiated as a function of a for infall from 
TojIM = 40, and as a function of Tq for six selected values of a. We also tabulate the total 
energy radiated (the sum of the I = 2, 3 and 4 modes) for selected values of Tq and a. The 
radiated energy was calculated using two different methods: direct numerical integration of 
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Eq. ( |2.9|) , and numerical integration of Eq. ( |3.7|) . In this way, the accuracy of the algorithms 
used could be tested. The results were found to agree to better than 1%. For comparison, 
we include in our figures the energy radiated by a particle falling from infinity, as calculated 
by Davis, Ruffini, Press, and Price (DRPP)||5T|. The DRPP result for the energy radiated 
in the / = 2 mode is {2M/fi'^)E2 = 1.84 x lO^^^ 
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FIG. 11: The total energy radiated in the / = 2 mode as a function of a for a particle falling in 
from ro/2M = 40. Outside of the range < a < 7.5, the total energy radiated exceeds the DRPP 
result. The energy is minimized for a ~ 1. 



We have seen that for large initial separations, the different contributions to the gravi- 
tational radiation do not interfere. The energy emitted in each multipole is then the direct 
sum of the energy contained in the initial-data pulses, the energy emitted in the acceleration 



radiation, and the energy radiated in the quasi-normal modes. In Sec. |III B| we showed that 
all choices of initial data result in two pulses at early times: the outgoing pulse and the 
backscattered ingoing pulse. These pulses have the smallest amplitude when a = 1. Hence, 
for large initial separations, the waveforms obtained with a = 1 are the ones that carry the 
least amount of energy to infinity. 



The total energy radiated in the / = 2 mode is shown in Fig. [11] as a function of a, for 
ro/2M = 40. We see that the energy is minimized for a ~ 1, and that it exceeds the DRPP 
result if a is outside the interval ^ a ^ 7.5. For such a large value of the initial radius, we 
confirm the general belief that the radiated energy is minimized if the initial hypersurface 
is conformally fiat[|g, |g. 

In Fig. |12| we plot the total energy radiated in the / = 2 mode as a function of Tq for 
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FIG. 12: The total energy radiated by the / = 2 mode as a function of ro- Displayed are six curves 
corresponding to the values a = 0, 0.55, 0.9, 1, 1.1, and 2. For these small values of Tq, a = 1 no 
longer minimizes the energy radiated. Instead, the value of a that minimizes the energy changes 
with Tq. This reflects the interference between the initial data, the particle, and the quasi-normal 
ringing contributions to the radiation. We were unable to find a single value of a that minimizes 
the energy in the interval 1 < ro/2M < 5. 

selected values of a. In the range 1 < ro/2M < 5, the energy is no longer minimized by 
the choice a = 1. Instead, we find that in the interval 2.7 < ro/2M < 5, the energy is 
minimized when a = 0.55, while in the interval 1 < ro/2M < 2.7, the minimum is achieved 
when a = 0.35. (These values are approximate, as we found it difficult in practice to locate 
the true minimum in the energy for a given r^ and /.) We also find that for a given r^,, 
no value of a minimizes the energy radiated for all modes. For example, for infall from 
ro/2M = 5 there is a minimum in the energy for the / = 3 and 4 modes when a = 0.9, but 
this value of a does not minimize the energy in the / = 2 mode (see Table |). 

Although we are unable to find a single value of a that minimizes the energy radiated 
for all / and over the whole interval 1 < r^/^M < 5, we see that in the cases considered, the 
total energy radiated is never minimized by choosing initial data with a > 1. Instead, for 
infall from a small Tq, it is minimized by a choice of a < 1. This is due to strong interference 
effects for initial data with a < 1 when Tq in the strong-field region of the spacetime; this 
interference was discussed in Sees. IIIIB| and [IIIC 



In Table | we display the total energy radiated in the Z = 2, 3, and 4 modes of the 
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TABLE I: Total energy radiated (in units of IMj [}}) for the ^ = 2, 3, and 4 modes, and a = 0.75, 
0.8, 0.85, 0.9, and 1.0. The sum of the first three multipole moments is denoted by E. For infah 
from short and intermediate distances, the minimum in E is achieved for a < 1. 





ro/2M 


I 


Q=0.75 


a =0.8 


a=0.85 


a=0.9 


a=1.0 


10.0 


2 


1.467e-2 


1.464e-2 


1.463e-2 


1.464e-2 


1.469e-2 




3 


1.961e-3 


1.876e-3 


1.820e-3 


1.789e-3 


1.790e-3 




4 


4.835e-4 


3.778e-4 


3.045e-4 


2.595e-4 


2.406e-4 




E 


1.711e-2 


1.690e-2 


1.675e-2 


1.669e-2 


1.672e-2 


5.0 


2 


1.305e-2 


1.309e-2 


1.316e-2 


1.326e-2 


1.353e-2 




3 


1.725e-3 


1.610e-3 


1.546e-3 


1.529e-3 


1.606e-3 




4 


5.600e-4 


3.883e-4 


2.761e-4 


2.156e-4 


2.249e-4 




E 


1.534e-2 


1.509e-2 


1.498e-2 


1.500e-2 


1.536e-2 


3.0 


2 


1.221e-2 


1.246e-2 


1.274e-2 


1.305e-2 


1.374e-2 




3 


1.359e-3 


1.307e-3 


1.323e-3 


1.398e-3 


1.693e-3 




4 


5.285e-4 


3.298e-4 


2.147e-4 


1.721e-4 


2.681e-4 




E 


1.410e-2 


1.410e-2 


1.428e-2 


1.462e-2 


1.570e-2 


2.1 


2 


1.238e-2 


1.310e-2 


1.383e-2 


1.457e-2 


1.608e-2 




3 


1.270e-3 


1.325e-3 


1.446e-3 


1.625e-3 


2.120e-3 




4 


4.134e-4 


2.515e-4 


1.757e-4 


1.748e-4 


3.598e-4 




E 


1.406e-2 


1.468e-2 


1.545e-2 


1.637e-2 


1.856e-2 


1.1 


2 


9.158e-4 


9.068e-4 


9.066e-4 


9.144e-4 


1.099e-3 




3 


2.017e-4 


1.892e-4 


1.854e-4 


1.894e-4 


2.250e-4 




4 


2.286e-4 


5.415e-5 


4.654e-5 


4.529e-5 


5.966e-5 




E 


1.346e-3 


1.150e-3 


1.139e-3 


1.149e-3 


1.384e-3 



Zerilli-Moncrief function, for selected values of a and r^- In general, the features presented 
previously for the case / = 2 are also present for higher multipole moments: the total energy 
radiated in a given multipole is not minimized by the conformally-flat choice of initial data 
{a = 1), but by some a < 1. We find that in some cases, the interference effects become so 
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important that the energy radiated increases with the multipole order /, e.g. E2 < E^ < £"4. 
This increase in energy with increasing / contradicts our slow-motion expectation, according 
to which the energy radiated should decrease with increasing multipole order. A typical 
example of this phenomenon is displayed in Table | for ro/2M = 1.1 and a = 0.75: here we 
find Ei> E^. 

IV. CONCLUSION 

A one-parameter family of time-symmetric initial data for radial infall of a particle into a 
Schwarzschild black hole was introduced in perturbation theory. The family is parameterized 
by a, which measures the gravitational-wave content of the initial surface. The family 
contains a conformally-fiat initial three-geometry as a special member [a = 1). Varying a 
allowed us to examine the infiuence of the initial data on the gravitational waves emitted 
by the particle-black hole system. 

We showed that for large initial separations, three stages can be clearly identified in the 
radiative process: initial-data-produced pulses, particle-produced acceleration radiation, and 
black-hole-produced quasi-normal ringing. For smaller separations, the three stages become 
confused and interference takes place. 

For large initial separations we confirmed the general belief that a conformally-fiat initial 
three-geometry minimizes the radiated energy. But we showed that for ro/2M < 10, the 
initial configuration that minimizes the energy is not the conformally-fiat choice. Instead, 
the configuration that minimizes the gravitational- wave content has a < 1. 

Most importantly, our numerical simulations show to what extent the gravitational wave- 
forms are infiuenced by the choice of initial data. As long as the particle falls toward the 
black hole from a distance ro/2M > 10, the part of the waveform that is associated with 
acceleration radiation and quasi-normal ringing is insensitive to the choice of initial data. 
In these cases, waveforms obtained with conformally-fiat initial data are an accurate rep- 
resentation of the true radiative process in the region of interest, because the unphysical 
radiation coming from the initial data propagates to null-infinity before the physical radi- 
ation becomes important. For ro/2M ^ 10, however, the two contributions interfere, and 
uncertainties associated with the choice of initial data hopelessly contaminate the waveforms. 
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APPENDIX: CONVERGENCE OF THE NUMERICAL CODE 

We now describe the convergence properties of our numerical algorithm. The method 
described in Sec. p.llA| is designed to be second-order convergent, i.e. the numerical solution 



converges towards the exact solution quadratically. This means that for a given choice of 
To, a, and I, the field at a fixed spacetime point is a function of A described by 

^N{A)=ijA + A^P, (A.l) 

where tpN is the numerical field obtained with a stepsize A throughout the grid, tpA is the 
exact solution, and p is an error term independent of A. Note that error terms of order 
(9(A'^) are neglected in this equation. 

The convergence of our numerical code can be tested by defining 

64j{A) = ^jv(2A) - ^^^(A) = SA^p, (A.2) 

and evaluating Sip for various values of A. This function is proportional to the error made in 
discretizing Eq. ( |2.3| ) and it satisfies SiplnA) = n'^6ip{A). This later property is important as 
we use it to determine the convergence rate of our code. We construct Sip from ip calculated 
with A = 0.005, 0.01, 0.02, and 0.04 on the null line v/2M = 500, i.e. at each grid point 
on this ingoing null line we calculate 6tp{0.005), 6iIj{0.01), and 6tp{0.02). In this way, we 
construct 6iIj{A) as a function of retarded time u. Note that for u < — r*, we are testing 
the vacuum finite-difference algorithm, whereas for u > — r* we are testing the convergence 
properties of the code involving the cells crossed by the particle's world line. 

From the definition of 6ip, it is obvious that they should satisfy 6ip{0.02) = ASiplO.Ol) = 
16(5?/'(0.005) at every grid point if Eq. ( |A.1| ) holds and the convergence is quadratic. In 
Fig. m, we display (5V'(0.02), 45^(0.01), and 16Sip{0.005), for the case ro/2M = 2.1, a = 1, 
and 1 = 2. The fact that the three curves are close together signals quadratic convergence, 
with 4(5^(0.01) and 165-0(0.005) being the closest. This is expected since as A is decreased, 
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FIG. 13: Convergence properties of our numerical code. We display 5^(0.02), 45^(0.01), and 
16(5^(0.005), for ro/2M = 2.1, a = 1, and / = 2. The agreement between the curves indicates that 
the numerical code is converging quadratically towards the exact solution. This is characteristic of 
the convergence in other cases. 



the approximation made in Eq. ( [A.l|) becomes more accurate, and we get a better coincidence 
between the curves. The convergence rate observed in this figure is typical of the convergence 
obtained for other choices of Tq, a, and /. 

By varying the location of the event horizon, and the location of future-null infinity, we 
have verified that our results are independent of the actual location of the null boundaries. 
We also verified that our results for a = 1 are consistent with the results presented by Lousto 



and Price lEH . 
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